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Abstract 



We study an electron that interacts with phonons or other linear or non- 
linear excitations as it resonantly tunnels. The method we use is based on 
mapping a many-body problem in a large variational space exactly onto a one- 
body problem. The method is conceptually simpler than previous Green's 
function approaches, and allows the essentially exact numerical solution of 
much more general problems. We solve tunneling problems with transverse 
channels, multiple sites coupled to phonons, and multiple phonon degrees of 

freedom and excitations. 
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We consider a single electron tunneling through a resonant tunneling diode or a quantum 
dot, in the presence of interactions with phonons or other excitations. This interaction leads 
to phonon assisted resonant tunneling and affects the peak-to- valley current ratio, which 
is important in device applications. Most previous treatments of this problem use a Green's 
function approach, often involving Keldysh formalism . An exact solution was obtained 
only for the special case of a single site coupled to a single phonon mode in Id 
Experimentally important transverse degrees of freedom were treated only to leading order 
in perturbation theory 0. 

The approach we use is to map the many-body problem in a possibly large variational 
space exactly onto a one-body problem. The Schrodinger equation, for the most part in real 
space, is then solved for the one-body problem. In our view, this is a conceptually simpler 
approach. It also produces explicit solutions for broad classes of problems that have not 
been solved before. The solution is essentially exact, in that the size of the variational space 
can be systematically increased until the answer converges. This approach need not make 
explicit use of Green's functions. (If desired, the Green's function can be recovered from the 
wavefunction.) Other inelastic tunneling problems can be solved by the same method. 

One can solve essentially any problem where a single electron tunnels, and where many- 
body interactions are limited to a finite region of space. The electron-phonon coupling may 
be nonzero on many distinct sites, including several sites in the quantum well, sites in the 
barrier, and sites in the leads near the barrier. The electron may couple diagonally or off- 
diagonally to many types of phonons, and multiple quanta can be excited in the same or 
different modes. The electron-phonon coupling may be nonlinear, and the phonons may 
have nonlinear interactions among themselves. The electron can have transverse degrees 
of freedom. Arbitrary one-body interactions, including barriers and disorder, can also be 
included. The electron can interact in an arbitrary way, including spin-flip scattering, with 
a group of interacting "captive" electrons in the tunneling region, so long as the captive 
electrons cannot escape into the leads. The method can also be used at non-zero tempera- 
ture. In practice, the method may make significant demands on computer resources when 
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more than about 10,000 inelastic channels are included. 
We consider the Hamiltonian 



H = Hel + Hph + Hel-ph, (1) 

Hei = J2 ^o^Vj ~ 51 + /i.e.), (2) 

3 j,k 

(3) 

m 

Hel-ph = -Yl \,rnc]cj{a\^ + a^)- (4) 



The potential on site j can describe a tunnel barrier, disorder, or a bias voltage. The 
hopping amplitude tj^k can vary from site to site. Xj^m is the (diagonal) coupling of an 
electron on site j to an optical or acoustic phonon mode m. A site can represent a single 



atomic Wannier orbital or a larger region of space |T0|. 

The method works for complicated barrier structures and interactions described by the 
above Hamiltonian, and for more general Hamiltonians. To illustrate the method in a simple 
context, however, we first consider the case with a single phonon mode that couples only 
to the electron density on site 0. The many-body problem is first restricted to a variational 
subspace. For illustrative purposes, only states containing 0, 1, or 2 phonon quanta are 
retained in the example. (A workstation could easily handle thousands of states.) The 
many-body scattering problem in the variational subspace is then mapped exactly onto a 
1-body problem with many channels, as shown in Fig. (|l|a) and explained in the caption. 
At zero temperature, an electron incident from the left is an incoming plane wave on the 
lower left lead. It has an amplitude to exit on any of the six leads, corresponding to elastic 
and inelastic backscattering and transmission. 

We seek the solution of the Schrodinger equation Eipj = J2k Hj^k'ipk on the tight-binding 
lattice of Fig. (|l|a) with the known eigenvalue E = —2t cos(/co) + e, where t and e are the 
hopping amplitude and diagonal energy of the left lead, and is the incoming wavevector. 
(The term "lead" refers to the translation invariant part of the system.) The boundary 
conditions are that an incoming wave is allowed on only one lead. The scattering problem 
is straightforward to solve. One method is to "prune" any lead that has only an outgoing 



wave. This exact procedure removes the lead from the problem, while changing the site 
energy on the last retained site to a value that is in general complex and changes with 
E. For example, the Schrodinger equation on site is Eipo = hi + eoipo — to^iipi, where 
hi includes the off-diagonal matrix elements from site to sites other than 1. To prune 
the lower right lead in Fig. (^), note that there is a unique outgoing wave of energy E 
in the lead, ipj = Aexp{ikij) . (If there are no propagating modes, choose instead the 
decaying mode ipj = ^6xp(— gij).) Back propagate this solution (through tunnel barriers 
if necessary) to obtain ipo and ipi, using the Schrodinger equation on each site j to obtain 
ipj-i as a function of ipj and ipj+i- The new Schrodinger equation on site 0, with all the 
sites in the lead removed, is now EipQ = hi + ioipo, where Cq = cq — ato,i) and a = ipi/ipo is 
generally complex. The leads to be pruned in Fig. (|l|a) are marked with vertical lines. The 
system to be solved after pruning is shown in Fig. @)). 

The problem in Fig. (|I|b) is so simple that it is best solved by a recursive trick, which 
does not work in general. To motivate the general solution, consider the same problem, but 
where the electron at site interacts with two distinct phonon modes of different frequencies. 
Again for illustrative purposes, choose a variational space that allows up to 2 phonon quanta 
in either mode, for a total of 3 x 3 = 9 phonon states (see Fig. |I|c). 

The pruned problem is not solved as a standard eigensystem, since the eigenvalue E 
is known in advance. Considering the amplitude ipo to be known, the problem is then 
to solve a system of (complex) linear equations of the form Ax = b, where x and b are 
vectors, with b proportional to ipQ. For this toy problem A is an 8 x 8 matrix. Once the 
system Ax = b is solved, the Schrodinger equation on sites 0, -1, etc., is used to determine 
the wavefunction on the first two sites of the left lead, and thus the coefficients oi and 
a2 in ipj = aiexp{ikoj) + 02 exp(— i/coj)- The current J leaving through the pruned leads, 
corresponding to elastic or inelastic transmission or backscattering in particular channels, is 
obtained using Jj^k = 2 Im(?/'^ tkj ipj). The formula is applied for a retained site j and a 
neighboring pruned site k. Current is conserved exactly, globally and at each vertex. This 
equation can be used to calculate ordinary current or a generalized current between two 
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many-body states |TT . 



Another form of electron-phonon coupling modulates the hopping matrix element t rather 
than the on-site energy e, 

H'el-ph = - J2 l3Amc]ck{a^n, + «m) + h.C. . (5) 

j,k,m, 

This off-diagonal coupling represents the fact that when an atom is displaced to the right, 
the hopping amplitude t to the atom on its right increases, because it is closer. A system 
with both types of electron-phonon coupling, which can be solved by the same method, is 
shown in Fig. (|l]d). 

Figures (2-4) show essentially exact results for problems that to our knowledge have 
not been previously solved. Figure plots transmission for more than one site coupled 
to a single phonon mode, and considers off-diagonal electron-phonon coupling. Figure 
has coupling to many distinct phonon modes of different frequencies [|12|. Finally, Fig. 



considers transverse degrees of freedom with electron recoil, which was previously treated 
only to leading order in perturbation theory 

As a simple test case. Fig. shows transmission through a quantum dot where a single 
phonon mode couples to an electron on a site with electron-phonon coupling strength A. 
Hopping matrix elements are t^^i = to between site and sites ±1 and tk,i = t for other 
nearest neighbors. We model the weak coupling through a tunnel barrier by a reduced to 
in this paper, although we could have just as easily used sites with increased tj. The same 
phonon mode also modulates the hopping matrix element (see Eq. (|^)). A variational space 
with up to 8 phonon quanta gives results accurate to the width of the plot lines. The inset 
of Fig. (^ shows the transmission for the case where phonons couple only to the electron 
density, i.e. 7 = 0. Our result agrees with previous calculations 0. The one-phonon 
sideband that we calculate was first seen experimentally by Goldman et al. Note that 
the low-energy "elastic" peak and the one-phonon sideband (at u ^ — X^/uJo) are each 
composed of both elastic and inelastic transmission. This is because there is an amplitude 
for the electron in the left lead to couple to the first excited state of the displaced harmonic 



oscillator for the electron on site 0, and then to tunnel into the right lead annihilating the 
phonon excitation, leaving the (now undisplaced) harmonic oscillator in its ground state. 
The calculation describes both inelastic tunneling, where phonons are emitted, and polaron 
physics, where phonons are emitted and reabsorbed. Elastic and inelastic tunneling can be 
experimentally distinguished by electroluminescence measurements W3 . 



Figures (^-d) show transmission when both diagonal and off-diagonal electron-phonon 
coupling is present. 7{l,_r} modulates the hopping between site and the {left, right} lead. 
The only difference between figures (|]a-d) is the relative sign of the coupling constants. 
There are clearly dramatic interference effects, which occur generically when a phonon is 
coupled to more than one site. The transmission peaks are much wider in Fig. (Qa) than in 
Fig. (0b). This is attributed to the fact that the line width is proportional to tg^, where the 
effective coupling to ~ + li^)- The oscillator acquires a positive displacement (x) due to 
the A coupling. Figures (0c) and Fig. (||d) represent the same physical system reversed left 
to right. The elastic part of the transmission is identical for the two cases, while the total 
transmission is quite different. This behavior is consistent with the unitarity requirements for 
a time-reversal invariant system with inelastic channels. The total integrated transmission 
for (0c) is quite different from (0d), because the sum rule is strongly violated for 7-type 
electron-phonon coupling. 

Figure (0) shows the transmission through a single site quantum dot, which is coupled 
to 10 phonon modes with different frequencies. The electron-phonon coupling is diagonal, 
as in Eq. (^). We compute transmission using variational spaces that contain up to a 
total Nph = 1, ... ,4 phonon quanta distributed in any way among the 10 modes. Using 



subroutines for large sparse systems [|14|, the system with Nph = 4 consisting of Nst = 1000 
states (2000 channels, including transmission and reflection) requires 20 seconds of CPU 
time per energy point on a Sparc 10 workstation. Each frequency point requires the solution 
of a Nst X Nst sparse system of complex linear equations. The inset of Fig. (^ displays an 
enlarged portion of the transmission function in the region of predominantly one-phonon 
contribution. The weak electron-phonon coupling A = 0.5 is in the experimentally relevant 



regime, where the single phonon peaks dominate those due to multi-phonons. Even so, it 
is clear that one does not get an accurate description of inelastic tunneling in a variational 
space containing only single phonon excitations, and that Nph = 3 is required to achieve 
reasonable convergence. 

It is straightforward to include transverse degrees of freedom, for example to model the 
case where the tunnel barrier is an extended, perhaps planar structure. We investigated a 
model containing Ny parallel leads with periodic boundary conditions in the (transverse) 
y— direction. The Hamiltonian, written with real-space indices j, / in the x— direction and 
momentum space indices fc, q in the y— direction is 

j,k j,l,k 

+ Y^kala,, - A/ cikCoMqH + ^-5)' (6) 

k k,q 

where the on-site energies are ej = [ei, eo, e^] for [j < 0, j = 0, j > 0] respectively. Hopping 
matrix elements are tj^i = for nearest neighbor j, I 7^ 0, and tj^i = t^o when j = or / = 0. 
Similarly tyj = ty for j 7^ and tyj = tyo otherwise. Diagonal electron-phonon coupling 
is restricted to sites where j = 0. Due to the translational symmetry in the y— direction, 
the total transverse momentum k is conserved. The electron momentum changes only as a 
consequence of the electron-phonon interaction. 

Figure @ shows the transmission for Ny = Q parallel leads. Figure (H) uses a variational 
space with up to Nph = 5 phonon quanta in any phonon modes, which gives satisfactory 
accuracy in the whole frequency and transverse momentum range. For simplicity we use a 
dispersionless phonon spectrum; however, generalization to momentum-dependent phonon 
frequencies and electron-phonon coupling is straightforward. At small transverse momentum 
A; = and tt/3, a strong nearly elastic resonance is located just below the noninteracting 
resonance at uj{k) = Cq — 2tyQCOsk. Inelastic side-peaks at higher u correspond to phonon 
creation, usually accompanied by a change in the electron momentum. For large transverse 
momentum k = 27t/3 or vr, resonant states where the electron creates a phonon and recoils 
to lower momentum can have a lower energy than the state with no phonons. This results 
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in side-peaks at energies below the weakened central peak. 

Generalizations: Finite temperature problems can be solved by having the incoming 
electron arrive on different leads in Fig. (|l|a) with the appropriate Boltzmann weights. It 
should now be possible to model more realistic coupled electron-phonon systems, to address 
such questions as why the inelastic peak is observed experimentally at the barrier phonon 
frequency rather than that of the well [|l],|l5|. It would be interesting to study the stronger 
electron-phonon coupling that arises when electrons are localized on impurities, or from 
phonon modes caused by crystal defects. It should also be possible to model more compli- 
cated band structures that include several Wannier functions per unit cell, and amorphous 
barriers if a suitable tight-binding description is known. 
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FIGURES 

FIG. 1. (a) Each dot represents a basis state \j,n) in the many-body Hilbert space. The 
rows of dots are the sites j = —3, . . . , 3 with n = 0, 1, or 2 phonon quanta. The bonds represent 
non-zero off-diagonal matrix elements in the Hamiltonian. The horizontal bonds are the hopping 
amplitudes tj^k- The vertical bonds represent the electron-phonon interaction, with amplitudes —A 
on the lower and — A\/2 on the upper vertical bond. The dots can also be interpreted as Wannier 
orbitals in an equivalent 1-body tight-binding model, (b) The system that results after pruning 
all but the incoming (lower left) lead. Sites with complex diagonal energies €j are shaded gray, 
(c) The pruned version of the problem where an electron at site couples to two distinct phonon 
modes. In the 3x3 grid of gray dots, the vertical direction is the number of quanta ni = 0, 1, 2 
of type 1 phonons, and the direction into the page is the number of quanta n2 of type 2 phonons. 
Before pruning, there was one right and one left lead attached to each gray dot. (d) A system 
with a single phonon mode, diagonal electron-phonon coupling of the type in Eq. on site 0, 
and off-diagonal electron-phonon coupling of the type in Eq. (|5|) between site and sites ±1. The 
pruned version is shown. 

FIG. 2. Transmission probability as a function of the incident electron energy. The heavy line 
is the total transmission, and the lighter (lower) line the elastic part. A single phonon mode is 
coupled to the electron density at site with strength A, and to the electron hopping amplitude 
between site and sites ±1 with strength 7. The other parameters of the Hamiltonian are: 
= Er = Co = 0.0 (no voltage drop across the dot), to = 0.2, luq = 1.0 and Nph = 8. All energies 
are in units of the hopping t in the leads. 

FIG. 3. Transmission through a single site that is coupled to ten phonon modes with phonon 
frequencies ujn uniformly distributed between uJmin and uJmax, with electron-phonon coupling con- 
stants A„ = 0.5u;„,. Different curves represent runs with basis sets containing a different total 
number of phonon quanta Nph- 
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FIG. 4. Transmission through Ny = 6 transverse channels as a function of the energy lo of the 
incident electron at four different choices of transverse momentum k, with txo = 0.2. The upper 
and lower lines represent total and elastic transmission respectively. 
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